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Abstract 

In [TH] the authors developed a method for computing normal forms of dynamical systems 
with a coupled cell network structure. We now apply this theory to one-parameter fami- 
lies of homogeneous feed-forward chains with 2-dimensional cells. Our main result is that 
Hopf bifurcations in such families generically generate branches of periodic solutions with 
amplitudes growing like ~ |A| 5 ,~ A| 5 ,~ l-M 15 , etc. Such amplified Hopf branches were 
previously found in a subclass of feed-forward networks with three cells, first under a normal 
form assumption |15| and later by explicit computations [8], [13] • We explain here how these 
bifurcations arise generically in a broader class of feed-forward chains of arbitrary length. 



1 Introduction 



In this paper we shall consider systems of ordinary differential equations of the following 
homogeneous feed-forward type: 
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f(x„ , x„-i, ■ ■ ■ , xi, xo; A). 



Here n + 1 £ N is the length of the feed- forward chain, the state variables xq, . . . , x„ 
are elements of a finite dimensional vector space V and the function / : V n+1 x R — s 



e v 

V is 

a parameter dependent response function. We shall assume that /(0; A) = for all A, and 
hence that equations (|l.ip admit a fully synchronous steady state solution x — (0, . . . , 0) for 
all values of the parameter. 

We are interested in the periodic solutions that emanate from this synchronous steady 
state as A varies. In order to find such synchrony breaking bifurcations of periodic solutions, 
let us denote by ai = Oi(A) := Di/(0; A) : V — I V. Then the linearization matrix of (|1.1[) at 
the synchronous equilibrium reads 
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This matrix displays a remarkable degeneracy: the eigenvalues of the submatrix oo each 
occur at least n times as eigenvalues of the linearization matrix (|1.2[) . Although such a 
1 : 1 : ■ • ■ : 1 resonance is highly uncommon in differential equations without any special 
structure, it is generic in the context of our feed-forward network. 

Assume for example that dim 1 !/ = 2 and that the eigenvalues of oo(A) form a complex 
conjugate pair that crosses the imaginary axis at nonzero speed. Then one may expect a 
very particular n-fold Hopf bifurcation to take place in equations (|l.ip . Theorem 11.11 below 
describes this bifurcation. It is the main result of this paper. 
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Theorem 1.1 Let V — R 2 = C and f : V n+1 x R -> V a smooth function. Assume that 
/(0; A) = and that Oo(0) = Do/(0; 0) : V — > V has eigenvalues Uluo 0. Then under 
otherwise generic conditions on f(X;X), the local normal form of fl.l]) near (x;X) — (0; 0) 
supports a branch of hyperbolic periodic solutions of the form 

x (t) = 0, xi(t) = fl 1 (A)e iw(A) * J . . . ,x n (t) = B n (X)e^ Wt 

of freguency uj(X) — ujq + 0(\\\) and amplitudes |Bi(A)| ~ |A| Ki , where Ki := ^ . 

The peculiar Hopf branch described in Theorem 1 1 . 1 1 has an xi-component with amplitude 
growing at the rate ~ |A| 5 of the ordinary Hopf bifurcation. The amplitude growth of its 
X2, X3, .. .-components at rates ~ |A| B ,~ | A| t» , . . . is much faster though. One could say 
that the feed-forward chain acts as an amplifier and it has been conjectured that this is why 
feed- forward chains occur so often as motifs in larger networks [14] . 

Due to their hyperbolicity, the periodic solutions of Theorem 11.11 persist as true solutions 
of (II. ip . We also remark that the branch of periodic solutions given in Theorem 11.11 implies 
that there exist many more branches of periodic solutions. This is because equations (|1.1[) 
admit a symmetry: whenever t H> (xo(t), . . . , x n (t)) is a solution to (|l.ip . then so is 

t H> (x (t),X (t),Xl(t), . . . ,X n -l(t)) . 

As a result, the Hopf bifurcation of Theorem 11.11 generates a total of n branches of periodic 
solutions (counted excluding the fully synchronous branch x(t) = 0), given by 

x (t) = . . . = xr-i(t) = 0, xr(t) = B 1 (\)e^ wt , . . .,x n {t) = B n ^ r+1 (\y^ x)t . 

Each of these branches emanates from the bifurcation point (x; A) = (0; 0) and they all have 
a different synchrony type. We will see that only the branch described explicitly in Theorem 
11.11 (the one with the largest amplitude and the least synchrony) can possibly be stable. 

Dynamical systems with a coupled cell network structure have attracted much attention 
in recent years, most notably in the work of Field, Golubitsky and Stewart and coworkers. 
For a general overview of this theory, we refer to [9], [12], [15], [20] and references therein. It 
has been noticed by many people that networks may robustly exhibit nontrivial synchronized 
dynamics [2], [3], [S], [S], |21j . |23j . |24| . Synchronous solutions may moreover undergo 
bifurcations with quite unusual features. Such synchrony breaking bifurcations have for 
example been studied in [I], [3], [B], [7], [TT] and [22]. 

The anomalous Hopf bifurcation of Theorem 11.11 has also been described in the literature 
on coupled cell networks [8], [13], [15], namely in the context of equations of the form 

x = g(x , x ; A) , 

xi = g(xi, x ; A) , (1.3) 
±2 = g(x 2 ,x 1 ;\). 

Note that (|1.3p arises as a special case of (|l.ll) . with n — 2 and / : V 3 — > V of the form 

f(X ,X u X a ) := 9 (X„,*i). 

In fact, the equivalent of Theorem 1 1.1 1 for equations of the form (|1.3p was first proved in [IS] 
under a normal form assumption. More precisely, it is assumed in [15] that the right hand 
side of (|1.3p commutes with the symmetry (3:0,^1,2:2) i-> (xo, e ls x\, e %a X2) and hence that 
g(X;X) has the symmetry 

g{e is X ,e is X 1 ;\) = e is g{X , Xi ; A) for s G R. (1.4) 

With this invariance, the periodic solutions of (|1.3|1 become relative equilibria, which makes 
them computable. We remark that a generic g(X;X) of course does not satisfy (|1.4p , but 
the existence of the anomalous Hopf branch was later confirmed for a generic g(X;X) by 
an explicit computation of the normal form of (j 1 .3 P in [8]. Finally, with the use of center 
manifold reduction, an alternative and more efficient proof of the same result was given in 
[13] . For similar normal form computations of other networks, we refer to [10] . 

Theorem 11.11 is thus an extension of the results in [8] , [13] and [15] , valid for a more 
general class of feed-forward chains of arbitrary length. More importantly, our proof of 
Theorem 11.11 is in the spirit of a generic bifurcation theory for coupled cell networks. It 
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makes use of the theory developed by the authors in [15] that explains how to compute the 
local normal form of a coupled cell network. 

In fact, we show in [18] that any network architecture admits a natural Lie algebra that 
captures the structure and invariants of the dynamical systems with this architecture. This 
structure is that of a so-called "semigroup network" and the main result in 18 is that the 
normal form of a semigroup network is a semigroup network as well. It turns out that 
equations (|1.1|) form an example of a semigroup network, and hence the normal form of 
(|1.1[) near a steady state is again of the form (|l.l|l . This proves it justified to assume that 
equations (|1.1|) are in normal form, and hence that f(X; A) satisfies 

f(e ls X , e ls X n ^,0; A) = e"f(X , X n -x,0; A) . (1.5) 

Using the latter invariance, the analysis of the periodic solutions of fll.ip is similar to the 
analysis in [15]. This analysis eventually leads to Theorem 11.11 It will also become clear 
that it is not reasonable to assume that the normal form of (|1.3|l is of the form (|1.3|l , because 
equations (|1.3|) do not form a semigroup network. 

This paper is organized as follows. In Section [2] we illustrate the amplifying character 
of our feed-forward chain by studying a synchrony breaking steady state bifurcation in case 
dim V = 1. This helps to fix ideas, and it yields an extension of some results obtained in [18j . 
In Section [3] we recall some results from [18] on semigroup coupled cell networks. Sections [4] 
and [5] are devoted to the algebra of linear feed- forward networks. The goal of these sections 
is to obtain, by means of linear normal form theory, a decomposition in semisimple and 
nilpotent part of the linearization matrix (|1.2[) . This is an essential preparatory step for the 
computation of a nonlinear normal form. We provide a linear normal form theorem in Section 
[3] and prove it in Section [5] using the concept of a semigroup ring. Finally, in Section H3 we 
use the .SW-decomposition obtained in Sections [4] and [5] to compute the nonlinear normal 
form of (|1.1[) under the assumptions of Theorem 11.11 A singularity analysis then leads to 
the proof of this theorem. 



2 An amplified steady state bifurcation 

In this section we describe a synchrony breaking steady state bifurcation in the feed-forward 
network (|1.1[) that confirms its amplifying character. This section is meant as an illustration 
before we start the actual proof of Theorem ll.il Most notably, we do not make use of normal 
forms in this section yet. 

We first of all remark that when /(O; 0) = 0, then equations (|1.1[) admit a fully syn- 
chronous steady state x = (0, . . . , 0) at the parameter value A = 0. This steady state 
persists to a nearby synchronous steady state (s(A), . . . , s(A)) for A close to under the 
condition that 

— /(s, . . . , s; 0) = ao(0) + . . . + o n (0) is invertible. 

as „ 

s— 

Throughout this paper, we will assume that this condition is met, so that by translating to 
the synchronous equilibrium we justify our assumption in the introduction that /(0; A) = 0. 

We are interested in branches of solutions that emanate from (x; A) = (0; 0). These are 
defined as follows: 

Definition 2.1 A branch of steady states near (a;; A) = (0; 0) is the germ of a continuous 
map x : [0, A ) — > V n+1 or x : (— A , 0] — > V n+1 with the properties that :r(0) = and such 
that each a;(A) is a steady state of the differential equations at the parameter value A. A 

To characterize the growth of a branch of steady states, we introduce a special symbol: 

Definition 2.2 For a branch x(X) of steady states we write a;(A) ~ |A| K , with k > 0, if 
there is a smooth function x* : [0, |Ao| K ) — ► V n+1 such that 

x(X) = x*(\\\ K ) and lim t^O- 

|A|no |A| K 

A 

The following theorem describes the branches of steady states that can bifurcate in (II. 1[) in 
the case that dim V = 1. It is a generalization of results in [18| . where it was assumed that 
n — 2. We spell out the proof of Theorem 12.31 in great detail and hope that this makes the 
structure of the proof of Theorem 11.11 more understandable. 
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Theorem 2.3 Let V = R and f : V" +1 x R — >• V a smooth function. Assume that /(0; A) — 

and i/iat Oo(0) = 0. Furthermore, assume the following generic conditions on f(X; A); 

1. o (0) + ... + o n (0) #0. 

2. ^(0)^0. 
5. oi (0) # 0. 

4. gf(0;0)^0. 

Then equations Sl.lp support In branches of hyperbolic steady states ( counted excluding the 
trivial steady state x = 0) near (x; X) = (0;0). More precisely, if we define m := jT=t f or 

1 = 1, . . . ,n, then for each r — 1, . . . , n there are two steady state branches with asymptotics 

X = . . . = X r -\ = 0, X r ~ |A| K1 , . . . , X n ~ |A| Kn - r+1 . 

Proof: Let us write Oq(0) = ^(0). We will give the proof of the theorem under the assump- 
tion that od(0)oi(0) > 0. The proof is similar when od(0)oi(0) < but it is inconvenient 
to give both proofs simultaneously. This technical problem will not occur in the proof of 
Theorem O 

Because Oo(0) + . . . + a n (0) ^ we can assume that /(0; A) = 0. Together with the 
remaining assumptions of the theorem, this implies that we can Taylor expand 

f{X; A) =a' o (0)XXo + ai(0)Xi + . . . + a n (0)X n + CXl 

+ O (|X„| 3 + |X | 2 |A| + |X ||A| 2 + ||(Xi, . . . ,X n )|| 2 (2.6) 
+ ||(X 1 ,...,X n )||-|A| + |X () |-||(Xi,...,X n )||) , 

in which 0o(0) + 0, oi(0) ^ 0, Oi(0) + . . . + o„(0) + and C := |^4(0; 0) / 0. With this in 
mind, we will solve for steady states by consecutively solving for xj the equations 

f(xj,Xj-i, . ..,xi,xo, ...,x ; A) = for j = 0, ...,n. 

First of all, since 0i(0) + . . . + O„(0) / 0, it holds by IpTo) that 

f(x ,. ..,x ; A) = (oi(0) + . . . + a„(0))a;o + O(|a;o| 2 + |A||x |) . 

We conclude that xq = := must hold for a steady state. 

In fact, it is clear for any 1 < r < n that xo = Xq = 0, . . . ,x r -i = x*_ 1 := provide 
solutions to the equations f(xo, ■ ■ ■ , Xo; A) = 0,...,f(x, — i,...,xo;A) = 0. Given these 
trivial solutions, let us try to find the possible steady state values of x r by solving 

/(av,av_i, ...,Xq;X) = ad(0)Ax r + Cxi + 0(\x r f + |x r | 2 |A| + |x r | ■ |A| 2 ) = . (2.7) 

Not surprisingly, one solution to this equation is given by x r = 0, but let us ignore this 
solution. Thus, we divide equation (|2.7|l by x r and consider the case that x r solves 

ad(0)A + Cx r + 0{\x r \ 2 + \x r \ ■ |A| + |A| 2 ) = . 

Because C ^ 0, the implicit function theorem guarantees a unique solution x r to this latter 
equation, given by 

^ = „^) A + 0(A 2 ). 

Note that x r is defined for both positive and negative values of A. Thus, we find two 
functions xf* : [0, A ) -> V such that x r {\) = xf*(\\\) ~ [A| i = |A| K1 solve (|2.7p . In 
particular, taking r = n, we proved the existence of two of the branches in the statement of 
the theorem, namely those with asymptotics xo = . . . = x n -i = and x„ ~ |A| K1 . 

In case 1 < r < n — 1, we shall ignore the branch k^T*(|A|) defined for negative A: it leads 
to an equation for x r+ i that can not be solved. This follows from computations similar to 
those given below for the positive branch. We will not provide these considerations in detail. 
We shall write x* = x^* ■ 

Concerning this positive branch, let us remark that linu^o = — an d have 

opposite sign due to our assumption that od(0)oi(0) > 0. This leads, for r < j < n — 1, to 
the following two induction hypotheses: 
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1. Assume that we found for all i = r, . . . , j certain smooth functions x* : [0, Ag 1 r+1 ) — > E 
so that 

xo=Xg = 0, . . . , a; r -i ~x*_ 1 = 0, x r = x*(X hl ) ~ A K1 , . . . , Xj = x*(X Kj ~ r+1 ) ~ \ n i-<-+i 
solve the equations /(xq, . . . , xq; A) = 0, . . . , f(xj,Xj-i, . . . , Xq, . . . , Xo; A) = 0. 

2. The sign of lim^o Xj ^K j _ T+1 - is opposite to the sign of . 

We remark that we just proved that these induction hypotheses are true for j = r. We will 

now try obtain Xj+i by solving the equation f(xj+i, Xj, . . . , Xq, . . . ,Xq; A) =0. 

Anticipating that the steady state value of Xj+i will be smoothly depending on X K i~ r + 2 , 

i 

let us at this point define the rescaled parameter fi := X Kj - r + 2 — and the rescaled 

unknown Xj+i =: fiyj+i. Then it holds t 
also the rescaled functions y* , . . . , y * by 



unknown Xj+i =: fiyj+i. Then it holds that ji = A which inspires us to define 



x*(X K1 ) =: (i y*([i ), . . . ,Xj (A Kj r+1 ) =: fj, y*j(fi ) for /i = A Kj-7 " +2 . 

By the first induction hypothesis, the functions y* are smooth and y* (0) = lim^o J xKj _ r+1 — 



7^ 0. Moreover, using (|2.6|) . one checks that in terms of the rescaled variables, the equation 
for Xj+i — nyj+i takes the form 

M 2 (a x (0)2/* (0) + Cy 2 +1 ) +OQy j+1 \ ■ | M | 3 + | M | 4 ) = 0. 

Dividing this by /i 2 , we find that we need to solve an equation of the form 

h(y J+1 ;v) = (ai(0)y*(0) + Cy 2 +1 ) + 0(\y j+x \ ■ \p\ + \fi\ 2 ) = 0. 

Recall that ai (0), y*j (0), C(0) 7^ 0. In fact, by the second induction hypothesis it holds that 

— ni ( '"^ J — > 0. Let Yj+x :— \J — ni ^°^> ^ so that ±Y,-+i are the solutions to the equation 
h(Y j+1 ,0) = ai(0)y*(0) + CY? +1 = 0. Then clearly D Vj+1 h{±Y j+1 , 0) = ±2CY j+1 ± and 
thus by the implicit function theorem there exist smooth functions yf^ip) = ±Y)+i + O(fi) 
solving h(yfj^ 1 (/j,) , ji) = 0. Correspondingly, the expressions 

x j+1 = (x* +1 ) ± (A K -'-+ 2 ) = mfti(n) = ±A Kj ~ r+2 Yj+i + 0(\ Kj ~ r+1 ) ~ A^-"+ 2 
provide two branches of solutions to the equation f(xj+i,x* . . . , x^; X) — 0. It holds that 

lim^o J+ ^ Kj _ T , +2 = J/j + i(0) = iiYj+i so that for precisely one of these branches the sign 

of this limit is opposite to the sign of . Only this branch can be used in the next step 
of the induction. This step is necessary precisely when j + 1 < n — 1. This finishes the 
induction and the proof that also for 1 < r < n — 1 there exist two steady state branches 
with asymptotics 

Xq — 0, ... , X r -1 — 0, X r ~ A M , . . . , X n ~ A*" . 

Note that these two branches only differ in their values of i„. 

We remark that in case a (0)ai(0) < 0, the branches exist for negative values of A. 

Finally, we consider the linearization matrix around a steady state on the r-th branch. 
This matrix is of the lower triangular form 

( Q (A) + ai(A) + ... + a n (A) ... \ 

* 00(A) ... 

* : 
00(A) 

: : * bi(A) 

\ * * ■■■ * bn-r+l W ) 



The eigenvalues of this linearization matrix are do (A) + . . . + a„(A), do (A) (r — 1 times), 
bi(A), . . . , b„_ r (A) and b n _ r+1 (A), where 

b 3 (A) = D f(x r+j _ u ...,xl) = 2Cx* +j _ 1 + OdAI"*- 1 ) ~ |A| K ^ . 
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for j = 1, . . . , n — r and 

bn-r+iW = D f(xt*, ...,x* ) = 2Cxt* + 0(|A|"»- ") ~ \\\^-r + i . 
For A / 0, these eigenvalues are real and nonzero, i.e. the branches are hyperbolic. □ 

Remark 2.4 Only one of the branches of solutions given in Theorem 12.31 can possibly be 
stable and this is one of the two branches with the least synchrony and the largest amplitude, 
i.e. one of the two branches with asymptotics 

X = 0, X\ ~ |A| K1 , . . . , X„ ~ |A| K ™ . 

This happens precisely when ao(0) + . . . + o„(0) < and ai(0) > 0. We leave the proof of 
this claim to the interested reader. A 

3 A semigroup network 

The feed-forward differential equations (|3.8p form an example of a so-called semigroup net- 
work. These networks were defined by the authors in [18], and they have the remarkable 
property that the Lie bracket of two semigroup networks is again a semigroup network. 

In the context of our feed-forward chain, this is perhaps best explained as follows. First 
of all, note that we can write the differential equations as 

x = Jf(x) for x £ V n+ , 
where the vector field 7/ : V n+1 -¥ V n+1 is denned for a function / : V n+1 — >■ V as 

(7f)j( x ) ~ f( x <T (j),x ai y), ■ ■ ■,%* n -- L {i),Xv n {i)) for j = 0, ...,n. (3.8) 
Here, uo, . . . , <x n are maps from {0, . . . , n} to {0, . . . , n}, given by 

di{j) := max{j -i,0}. 
One can now observe that for all < ii, 12 < n it holds that 

0~i±0~i 2 '•— 0~i± 0~i 2 = 0~ Tn in{i-L+i2 ,n} • 

This means in particular that the collection 

E := {(To, <Tl, . . . , <r n } 

is closed under composition: it is a semigroup. In fact, E is commutative and generated by 
the elements <to and o\. Moreover, the elements ctq and cr n are somewhat special: 00 is the 
unit of E and o n plays the role of "zero", because o n Oi = Uicr n = a n for all < i < n. 

All of this leads to the following result, that was proved for general semigroup networks 
in [18]: 

Theorem 3.1 Define for < i < n the linear map A <Ti : V n+1 -»■ V n+1 by 

(A tTi X)j := X m ^ n {j +i n y . 

Then it holds that 

A ail o Aa i2 = Aa H a i2 and that A ai (x ao (j-), . . . , x an{j) ) = (ajffoC^o)), • • • , ^(^(j))) ■ 
Moreover, for any f,g G C 00 (V™ +1 , V) it holds that 

hf,Jg] = 7[/,s] s 

where 

n 

[f,gh(X) := J2 °if( x ) ■ 9(A ai X) - D i9 (X) ■ f(A at X) . 

i=0 
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Proof: It is clear that 



(A tTii A tTi2 X)j — X miB ^j +il+i2in y — {A arainiii+i2 n yX)j — (Aa^a^X)] . 

The fact that A tTi {x tTQ ^j- ) , . . . ,x an ts)) = i x a (a (j))i ■ ■ • 1 x cnWU))) is obvious from our defini- 
tions. This proves the first claim of the theorem. 

Next, recall that [7/ ,7 9 ](a0 = D~/f(x) ■ J g (x) — D~/ g (x) ■ 7/(2;). One computes that 

n 

{D la {x) ■ lf {x))i =Y t D„{rtf(x)) i ■ (i a ) k (x) = 

n , 

fe=0 fc 
n 

y^gi/(x CT0 (j), . . .,x an(j) ) ■ g(x ao(rTiU)) , . . . ,s CT „( CTs (j))) = 

i=0 
n 

Dif(x ao(j) , . . . ,x anU) ) ■ g(A ai (x^ oU) ,.. -,x an( j))) . 

i=0 

With a similar computation for (D r y g (x) ■ "ff{x))j, we thus find that 

n 

hf,l g h(x)=^D i f(X).g(A ai X)-D i g(X)-f(A ai X) 
This proves the theorem. □ 



The first statement of Theorem 13.11 is that the map 

o-i^Aat, E -> gl(V n+1 ) 

is a representation of E. The second and more important statement is that the Lie bracket 
of the two feed-forward vector fields 7/ and 7 9 is another feed-forward vector field of the 
same form, namely 7[/, s ] E . Moreover, the new response function [/, <?]s is computed from / 
and g with the use of the representation o~i h> A ai . 

The most important consequence of Theorem 13. II is that the collection 

7 (C 00 (V m+1 , V)) = { 7/ I / G C°°(V n+ \ V)} 

of feed-forward vector fields is a Lie algebra. This implies for example that there exists a 
large class of transformations of the phase space V n+1 that leaves the class of feed- forward 
vector fields invariant: the time-1 flow e 7s of any feed-forward vector field 7 9 , will transform 
the feed- forward vector field 7/ into another feed-forward vector field, namely: 

(e 7 »), 7/ = e ad ^ ( 7/ ) = 7/ + [j g ,jf] + Ills, ha, If]] + ■■■ = 7/+[ 9 ,/ ]s +i [ s ,[ 9 ,/] s ] s+ ... • 

This explains why transformations of V n+1 of the form e 7s play an important role in the 
theory of local normal forms of semigroup networks. Theorem 13.21 below, for example, was 
proved in [15]. It describes the normal forms of one-parameter families of feed- forward 
networks. To formulate it, we define for k, I > 0, 

P k l : = {/ : V n+1 xR — > V homogeneous polynomial of degree k+1 in X and degree I in A} . 



Theorem 3.2 (Normal form theorem) Let f G C°°(V n+1 x R, V) and assume that 
/(0; A) = 0. Let us write the Taylor expansion of f as 

f = (/o,o + /o,i + /o.a + •••) + (A,o + + /i.a + ...) + ••• 
A := D x7f (0; 0) = 7/ ,o = A s + A N : V n+1 -> V n+1 



with f kll G P h ' 1 . We moreover denote by 
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the SN -decomposition of the linearization /(0;0) . Finally, let 1 < n,r2 < oo. 

Then there exists a polynomial family A n> $(•; A) of analytic diffeomorphisms, defined 
for A in an open neighborhood of and each sending an open neighborhood of in V n+1 to 
an open neighborhood of in V n+1 , such that $(•; A) conjugates 7/(.;A) £° Tf(-\)> where 

1 = (/o,o + 7 ,i + 7o, 2 + ■•■) + (7i,o + 7 M + 7i. 2 + •••) + ••■ 

/las t/ie property that for all < k < ri and < I < r2, 

e aAs °lj kl =lj kl ° e As for allseR. (3.9) 

Proof: [Sketch] This proof is based on the fact that the spaces P k ' 1 are graded, that is 

jpfe.i pif,ij s ^ pk+K,l+L 

As a consequence, one can start by choosing an go,i £ P ' 1 and use the time-1 flow e 790 1 of 
j go l to transform 7/ = 7/o,o+/o,i+ - into 

( e7so,1 )*7/ = 7/+[bo,i,/]e+- = 7/ ,o+(/o,i+[9o,i,/o,o]e)+- ■ 

By choosing go,i appropriately, one can then make sure that jj ^ :— jf 1 + [ go 1 j ] E com- 
mutes with the semisimple part As of jf g . This is a consequence of the fact that both As 
and An are semigroup networks. The precise argument leading to this result is nontrivial 
though and is given in detail in |181 Section 9]. 

Next, one normalizes / , 2 , . . . , /o,r 2 , /i,o, ■ • ■ , fi,r 2 , ■ ■ ■_ L fr 1 ,o, ■ ■ ■ , f ri ,r 2 ■ The grading of 
the P k ' 1 ensures that, once fk,i has been normalized into f kl , it is not changed anymore by 
any subsequent transformations. □ 

In short, Theorem 13.21 says that we can arrange that the local normal form of a parameter 
family of feed-forward vector fields 7/ is another parameter family of feed-forward vector 
fields 7j. Moreover, this normal form may be assumed to admit a "normal form symmetry" : 
it commutes with the flow of the semisimple part As of the linearization -0^7/ (0; 0) = 7/ • 



4 A linear normal form 



Theorem 13.21 describes the normalization of the parameter dependent coupled cell network 
vector field 7/ with respect to the semisimple part of the linearization 

A:= D xlf (0;0) = 7/0,0 • 

We recall that /o,o(^) = ao(0)Xo + ■ + a n (0)X n for certain 

o (0),...,a n (0) G :=0((V) 

and hence that the matrix of A is given by 



/ oo(0) + oi(0) 
01 (0) 



A = 



■o„(0) 
■o„(0) 





oo(0) 



V 



o„_i(0) + a„(0) a„_ 2 (0) 
On(0) a„_i(0) 



0o(0) 
01 (0) o (0) / 



(4.10) 



To determine the semisimple part of this matrix, we will bring A in "linear coupled cell 
network normal form" . This linear normal form is described in the following theorem: 

Theorem 4.1 Let A be the (n + 1) x (n+ 1) matrix given in Jj.Kfy and assume that ao(0) 
is semisimple. Then there exist linear maps gi, . . . ,g n -i '■ V n+ — >• V of the form 



gi(X) = bi(Xi - X n ) with bi e g 
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so that the consecutive time-1 flows of the linear maps 7 9i : V n+1 — > V n+1 

A := (e 7 ""-i o . . . oe 1j i) o Ao (e~ 7s i o . . . oe~ 7 ""-i) = 
/ oo(0) + ai(0) + ... + o„(0) •■• \ 



orm A into 



Oi(0) + 



0n(0) Oo (0) 



V 



On-l(0)+On(0) 5«-2(0) 

o„(0) a„-i(0) 



oo (0) 
fii (0) ao(0) / 



for which it holds that 

[oo(0),Bi(0)] = 0o(0)oi(0) - 0i(0)oo(0) = for all 1 < i < n 
As a consequence, A admits a decomposition A = As + An with 



As 



( oo(0) + Oi(0) + ... + q„(0) ••• 
Oi(0) + ... + o n (0) oo (0) ■•■ 



oi(0) + ... + o„(0) ■■■ oo(O) 
V Oi(0) + ... + On(0) ■■■ oo (0) j 



such that the map Am is nilpotent and 

[A s , A N ] = A s A N - A N A S = . 

The map As is semisimple if and only if oo (0) + 5i (0) + . . . + a n (0) is semisimple. 

We call the matrix A of Theorem 14.11 the linear almost normal form of the linearization 
matrix A. When Oo(0) and Oo(0) + 5i(0) + . . . + o n (0) are both semisimple, then the desired 
SW-decomposition of the linearization can be read off from this almost normal form. 

We shall prove Theorem 14.11 in Section [5] below. But before we do so, we would like to 
provide an alternative proof in case n = 2 here: in this case the theorem follows quite easily 
from an explicit matrix computation. The proof in Section [5] will be a bit more abstract. 

Proof (of Theorem 14.11 in case n = 2): We shall put A = and write Oi = Oi(0). Then, 
in case n = 2, the matrix (|4.10[) takes the form 



Oo + Oi + 02 








Ol + 2 


Oo 





Q2 


Ol 


Oo 



A = 



for certain oo, 01,02 € g. We can decompose A as A — As + An with 



00 + Ol + 02 
Ol + 02 

Ol + 02 




and A? 




(4.11) 



(4.12) 



It is clear that An is nilpotent. In addition, we can think of As as semisimple, because 

(idy, idy, idv), (0, idv, 0) and (0,0, idy) 

are "eigenvectors" of As with respectively the "eigenvalues" 

Oo + Oi + 02, ao and 00 . 

Of course, these eigenvalues are actually linear maps, namely elements of g, and it is clear 
that As is truly semisimple only when these eigenvalues are semisimple elements of g. 

The decomposition A = As + An may of course not be the SW-decomposition of A, 
because As and An in general need not commute. In fact, one computes that 



Oo + 01 + 02 

oi + 02 00 
oi + 02 Oo 




-[ao, oi 





[00, 01] 
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We shall resolve this problem by transforming A to a matrix A for which [do, Qi] = 0. This 
works as follows. First, we define, for some bi € g, the function 

gx : V 3 V by g 1 (X ,X 1 ,X 2 ) := bi(Xi - X 2 ) . 

Then j gi : V 3 — > V 3 is a linear map of the form 



791 




Moreover, it holds that 



exp 




A little computation now shows that conjugation with e 7si transforms A into 

_ / oo + m + a 2 \ / ao + ai + 02 \ 

A = e 7sl o A o e~ 7sl = j oi + o 2 Qo = Si + a 2 a . 

\ 02 + [ao,bi] ai-[a ,bi] a / \ a 2 Oi o / 

Here, we defined 

Si := ai - [o , bi] and a 2 := a 2 + [o , bi] . 

The essential step is now to choose bi in such a way that Qi commutes with ao. This is 
possible because we assumed that ao is semisimple in g, so that 

= im ad„ © ker ad ao in which ad 0o : b t-¥ [oo , b] , g — > o . 

Hence, we can decompose oi = oi m + Oi er for unique o™ € imad 0o and 0i cr € kerad 0o . If we 
now choose bi so that ad ao (bi) = oi m , then we obtain as a result that 6i = Cli — [do, bi] = 
oi — ad ao (bi) = oi — oi m = 0i cr , and hence that 

[oo, oi] = [oo, o^ cr ] = ad cl0 (o 1 ^ cr ) = . 

With such choice of bi, we obtain that As and Ajv commute as required. □ 

Remark 4.2 In the process of normalizing ai into 6i = Oi — [oo,bi], we automatically 
change a 2 into a 2 = a 2 + [ao, bi]. This means that o 2 will in general not be zero if a 2 = 0. 
Thus, already when we put the linear part of equations (|1.3[) in normal form, we obtain a 
system that is not of the restricted form (|1.3p but of the form A 



Our proof of Theorem l4.1l for general n is similar to the above matrix computation for n = 2. 
One could give this proof using the same matrix notation, but we found it more convenient 
to introduce a more abstract setting first. 

5 A semigroup ring 

In our proof of Theorem 14. II we shall make use of formal expressions of the form 

n 

aitJi with Qi G g and ai £ E . 

i=0 

We shall denote the space of such expressions as 



Moreover, when 



M = g£ = j^OiCTi j ai £ g,cr 4 e E j 



A = sxiOi and B = bjCj 

i=0 j=0 
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are two elements of M, then we define their product in M as 

n n n n 

i=0 j=0 i=0 3=0 



(5.13) 



The collection M is an example of a semigroup ring, cf. |17l Example 1.4]. It can be vieved 
as a module over the ring g with basis E, and also as a representation of E. For us it will 
suffice to think of M as an associative algebra that inherits its multiplication from g and E. 
One can remark that this algebra is graded: if we define for < k < n the collection 

M k := g{a k , . . . , a n } = ^ (ucrj Ck G g > C M , 



then M = M and for any < k, I < n it holds by (f5TT5|l that 

in{fc + Z,n} ■ 

Thus, each M k is an ideal in M and 

M = Mo D Mi D • • • D M„ D 

is a filtration. 

M also has the structure of a Lie algebra. Using that E is commutative, we find that 
the Lie bracket of two elements of M is given by the relatively simple expression 



[A, B] M := A- B- B ■ A = J2 XK E>j]<W{i +J >}- 

i=0 j=0 

The role of the semigroup ring M is explained in the following proposition: 
Proposition 5.1 The assignment 



( oo + Qi + . . . + a„ 

Ol + - • • + On O0 



\ 





V 



Cln-1 + On On-2 
On On-1 



Q0 
Ql Oo 



is a homomorphism of associative algebras from M to gl(V n+1 ). 
Proof: For < i < n, let us define the (n + 1) x (n + 1) matrix 



/ o» 
o t 



o t 
o, 













V o 







\ 





o / 



(5.14) 



With this definition, o^ Ni is the matrix of the linear network vector field 7 (diXi), where by 
tXiXi : V n+1 —¥ V we denote the map X diXi. 

Moreover, the assignment of the proposition is given by 



M 3 J2 Oi0"i ^ a i N i £ 0'(^ n+1 ) 
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It is easy to compute that 

(tttNiXbjNj) = (a i b j )N inill{i+j!n} . 
By distributivity of the matrix product, it thus follows that 

/ n \ / rt \ tl tl 

I £ cuNi) ]T bjNj = ]T £(o<J'i)^inta{i+-j,»} • 

\i=0 / \j=0 / i=0 j=0 

This product is homomorphic to the product in the semigroup ring (|5.13[) . □ 

It clearly has notational advantages to represent matrices of the form (|4.10p by elements 
of the semigroup ring, so this is what we will do in the remainder of this section. We also 
choose to perform the matrix computations that are necessary for the proof of Theorem 14.11 
inside the semigroup ring M and not inside gl(V n+1 ). 

But we stress that Proposition 15.11 proves it justified to think of JZILo <Xi ° i ^ M as the 
(n + 1) X (n + l)-matrix JZILo °»-W*> which in turn is the matrix of the map 



Proof (of Theorem 14. 1J) : Let us define for < i < n, 

Hi := Gi — a n £ Mi . 

Then it holds that /i n = and that 

Ms/Aj ~ (°» ~ a n){°~j ~ 0~rt) = °"mm{j + i,n} — °~ n = Mmin{t+j,n} • 

In particular, every \n is nilpotent in M. We also remark that the collection 

flO, ... , (J.n-1, °~n 

is a g-basis for the module M. 

Let us now assume that for some 1 < k < n — 1 there exist elements of the form 

Gi = bifii,...,G k -i = b k -!fik-i £ M 

such that 

n 

A:= (V 3 *- 1 •■■e Gl ) • A • (V Gl ■■■e- Gfc - 1 ) =^5,(7, 

»=o 

has the property that 5o = ao and ad ao (oi) = [oo, at] = for all 1 < i < k — 1. 
Then we pick £ arbitrary and define 

Gk '■= bfe/ife . 

Clearly, Gk is nilpotent in M and hence e Gfc is a finite series expansion and easy to compute. 
One finds that 

e° k ■ A ■ e~ Gk = e ad «fc {A) 



oo 1 



OO ^ Tl 



3 =0 J ' i=0 
k 

= ^2 0,0-j - [a , b k ]o- k mod Affc+i. 



As was explained in Section [2] we can now choose bk in such a way that Ofe := St — [ao, bk] 
commutes with ao, because Oo is semisimple. 
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This proves, by induction, that A — X^j=o a i°j can De normalized into 

71 

A = Oqito + OtiOi such that [ao, Si] = for 1 < i < n — 1 . 
Next, let us decompose A = aocro + X^"=i ^ i<Ti as ~ ^ s + An with 

/ Tl \ n — 1 n — 1 

As = aocro + I a» 1 a n and Ajv = ^ Si (<Xi - a n ) = ^ Oi/ij . 

\i=l / i = l i=l 

Then, because o~„/Xi = 0, it is clear that [As, Am]m = 0, that A^ is nilpotent, and that 

Asfn = aofii for i = 0, • • • , n — 1 and Ascr n = I ao + ^2 K» ■ 

This shows that no, ... , fi n -i and <t„, form a basis of "eigenvectors" of As with respectively 
the "eigenvalues" ao,-- - , ao and ao + J^"=i in particular, As is semisimple if these 
eigenvalues are semisimple in a. □ 

Remark 5.2 If one wishes, one can also bring ao in Jordan normal form. Indeed, if bo G g 
is an invertible map for which boaobg - ls m Jordan normal form with respect to some basis 
of V, then one can define Go = booo. We then observe that Go is invertible in M, with 
Gq 1 = bg Vo, and that 

n 

Go ■ A ■ Gq 1 = ^(boOibjY 1 ) m . 

This shows that one may assume that the matrix of ao is given in Jordan normal form. A 

Remark 5.3 For the proof of Theorem l4.ll the assumption that ao is semisimple is essential. 
If ao is not semisimple, then it can in general not be arranged that [ao, a>] = and hence it 
is not clear that [As, An]m = either. A 

Remark 5.4 We stress that for As to be truly semisimple, it is necessary that its eigen- 
values ao, ao + Si + . . . + On are semisimple in g. Nevertheless, it will become clear that this 
is not important for us. A 

Remark 5.5 One could also try to normalize a n by conjugating with an element of the 
form e G " , with G„ = b„(r n for some b„ G g. This indeed transforms a„ further, but we 
remark that this normalization is: 

1. Much more difficult to carry out than the normalization of the Oi with 1 < i < n — 1, 
because o n is not nilpotent but idempotent. 

2. Not necessary for the computation of the ^^-decomposition of A. Recall that we are 
interested in this SW-decomposition for the computation of a nonlinear normal form. 

For us, the equivalent of the Jordan normal form of A = X]"=o aiGi wm therefore be a 
matrix of the form A = aooo + X]™=i Sicr^ with [ao, Si] = for 1 < i < n — 1. Since we do 
not normalize a„, we call A an almost normal form of A. A 



6 Analysis of the nonlinear normal form 

In this section, we prove Theorem 1 1.1 1 of the introduction, that is formulated more precisely 
as Theorem 16.11 below. This theorem generalizes the results in [8], [14] and [15] to feed- 
forward chains of arbitrary length. We like to point out that, except for the linear algebra, 
various computational aspects of our proof are the same as those given in |15] . 

Before we formulate Theorem l6.ll we remark that when V = R 2 and ao(0) : V — > V has 
nonzero eigenvalues ±iwo, then we may assume that ao(0) is given in Jordan normal form 



ao (0) = 




As a consequence, it is then convenient to identify X — (X 1 , X 2 ) G R 2 with the complex 
number X 1 + iX 2 G C. This turns the linear map ao(0) : R 2 — > R 2 into the multiplication 
by iu)o from C to C and allows for Theorem 16. ll to be formulated conveniently as follows: 
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Theorem 6.1 Let V — R 2 = C and f : V n+1 x R -> V a smooth function. Assume that 
/(0; A) = and that oo(0) has eigenvalues zkiojQ 7^ 0. Furthermore, assume the following 
generic conditions on f(X; A); 

1. Persistence of the steady state: Oo(0) + . . . + o„(0) is invertible. 

2. Eigenvalue crossing: ^ j A _ tr do(A) 7^ 0. 

3. Nilpotency: troi(O) / 0. 

4. Nonlmearity: ^^ffl + 0. 

Then the local normal form of equations $1.1]) near (x; A) = (0; 0) supports n branches of 
periodic solutions (counted excluding the trivial steady state x{t) = 0). They are given for 
r — 1, . . . , n and for certain Bi(X), . . . , B n (\) 6 C by 

x (t) = . . . = xr-i(t) = 0,x r (t) = S 1 (A)e iwWt ) . . . ,x n {t) = B n _ r+1 (A)e^ (A)t . 

Frequency and amplitudes satisfy u)(\) — ujq + C(|A|) and |_Bi(A)| ~ |A| Ki with m := ^gi^r 
fori = 1, . . . ,n. These branches are hyperbolic if and only i/ao(0) + . . . + o n (0) is hyperbolic. 

Before we prove Theorem 16. II let us summarize the results of the previous sections as follows: 

Proposition 6.2 Under the assumptions of Theorem \6.1\ the vector field 7/ admits a local 
normal form fj near (x; A) = (0; 0) for which it holds that 

J{e tbJ ° s X , e <l " o %_i,0; A) = e lbJ ° 3 J{X , X n -i,0; A) for allsSR. (6.15) 
Proof: By Theorem 14.11 there exist linear maps gt : V n+1 — > V (1 < i < n — 1) of the form 

gi(X) = bi(Xi - X n ) so that e^-i 0...0 e 7 " : V n+1 -> V n+1 
transforms A = D x ~/f(0; 0) = 7/ into 

A := (e 79 "- 1 o . . . o e 7si ) o A o (e~ 7si o . . . o e -78 "- 1 ) in "almost normal form" . 
This means that A decomposes as A = As + An with An nilpotent, [As, An] = and 

As(x) = (0, oo(0)a;i, . . . , oo(0)a;„) + O(x ) = (0, iu xi, iu x n ) + O(x ) ■ 
In particular, we find that the subspace 

V n+1 := {x £ V n+1 I xo = 0} 

is invariant under both As and An- 

It is clear that the restriction of As to V " +1 is semisimple. Because An is nilpotent, 
this proves that A = As + An is the STV-decomposition of the restriction of A to Vq + . 

Because the transformation e 79 "- 1 o . . . o e 7fll is a composition of flows of semigroup 
network vector fields, it transforms 7/ into another coupled cell network 

( e T3n-l o . . . o e 7fll )*7 / = Jj. 

It obviously holds that 

^7 7 (0;0) = 7 /oi0 =3. 

By Theorem l3.2l it can now be arranged that the local normal form of 7j near (x; A) = (0; 0) 
is a network 7j that commutes with the flow of As, i.e. 

7j(e sAs a;) = e sAs (7j(a;)) for all s G R . 
Restricted to V$ +1 , this means that 

7(e iWo %-,...,e iwos a;i J 0,...,0;A) = e iu ° s J{x h . . . , xi, 0, . . . , 0; A) for all 1 < j <n. 
The latter is true if and only if (|6.15|l holds. □ 
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Remark 6.3 We remark that on the invariant subspace V n+1 the normal form symmetry 
of 7j reduces to the classical normal form symmetry 

(0, Xi,... ,x n ) >-> (0, e l "° s xi, . .., e %u ° s x n ) 

of the 1:1:...: 1-resonant harmonic oscillator. One could also try to normalize yj further 
with respect to the nilpotent operator An, see [19] . We will not exploit this further freedom 
though. A 

We are now ready to prove Theorem 16.11 using computations similar to those given in [15] : 

Proof (of Theorem 16 . 1[) : In this proof, we shall use the normal form symmetry (16.15[1 . In 
fact, applied to a monomial 

f(Xo, . . . , X n -i, 0; A) = Xq° ■ ■ • xf"/ • X~l ■ ■ ■ X^Lj 1 ■ X s 

equation (|6.15p yields the restriction that 2j=o (Pi ~7i) = !• Hence, the general polynomial 
or smooth f(X; A) that satisfies (16.15P must be of the form 

n-l 

J(X , Xn-i,0; A) = Y, F i( X ^l M =o,...,»-i i X ^ 

for certain polynomial or smooth complex-valued functions Fj of the parameter and of the 
complex invariants 

X k X l = {Xlxf + X 2 k Xf) + i(X 2 k Xl - Xixf) with < k, I < n - 1 . 

Moreover, the assumptions of the theorem imply that the first two of these functions Taylor 
expand as follows: 

F {. . . , X k X u . . . ■ A) = iivo + qA + C*|X | 2 

+ O (|X | 4 + |A| ■ |X | 2 + |A| 2 + IKXr,.. .,X n )\\ 2 ) , 
Fi(. . . , X k X h ...;X)=/3 + O(\\(X , X n )|| 2 + |A|) . 

Here, we defined C := | diRcXo)* ^ ^' ^ or wmcn ^ h°lds by assumption that ReC 7^ 0. 
Moreover, a — ai + iai € C and ft = /3i + ifii € C are those numbers for which 

aiA + 0(|A| 2 ) -a; -a2A + O(|A| 2 ) \ 
ao(A) -{iuo + a 2 \ + 0(|A| 2 ) Ql A + 0(|A| 2 ) ) and 



01(A) 



Pi+0(\\\) -/3 2 + 0(|A|) 
fr+0(\\\) Pi + 0(\X\) 



In particular it holds that 



1 d 
ai= 2dX 



troo(A) = i — trao(A) / and that 



/3 1 = itro 1 (0) = itra 1 (0)/0. 

Here, we have used that the trace of So (A) = 00(A) — [oo(0), bo (A)] is equal to that of do (A) 
and the trace of Si (A) = ai(A) — [oo(0), bi(A)] is equal to that of 01(A). 

Using this information, we will now try to find the relative equilibria of the normal form. 
To this end, we set xo(t) = Xg(t) := and solve for the periodic solutions of (|1.1[) by solving 
for j — 1, . . . , n the equations Xj = f(xj,Xj-i, . . . ,Xi, xq, . . . , xq; X) — consecutively. 

To find xi(t), we solve 

^ =/(xi,0,...,0;A) = F (|xi| 2 ,0; \)xi = 

(iw +aA + C|a;i| 2 +O(|a:i| 4 + |A| ■ \x \ 2 + |A| 2 )) an . (6.16) 
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Equation ()6.16|) is the normal form of the ordinary Hopf bifurcation. Its relative equilibria 
are found by making an Ansatz xi(t) = B\e llJ}t for B\ G C and cj G R. Ignoring the solution 
Bi = 0, this Ansatz leads to a complex equation in |-Bi| 2 and w: 

iw = iuo + aX + C|Bi| 2 + 0(|Bi| 4 + |A| ■ |Bi| 2 + |A| 2 ) . (6.17) 

The real part of this equation 

= aiA + ReC|Bi| 2 + 0(|Bi| 4 + |A| ■ |Bi| 2 + |A| 2 ) (6.18) 

can only be solved for small values of A and for small but nonnegative values of |-Bi| 2 if 

aiA/ReC < 0. 

This criterion determines the sign of A and thus whether the bifurcation is subcritical or 
supercritical. Thus, from now on we choose the sign of A so that aiA/ReC < 0. 

Because there is no restriction on the argument of B\, let us try to find B\ = |_Bi| > 
real. Anticipating that |_Bi| will be a smooth function of |A|5, we define fi := |A| 3 and 
=: fiZi. Then equation (|6.18p becomes a smooth equation in Z\ and /i of the form 

li 2 (|qi| - \ReC\Z 2 )+O{fi 4 ) = 0. 

Dividing this equation by fi 2 , we see that we need to solve an equation of the form 

h{Zi,n) = |qi| - \ReC\Z 2 +0(fi 2 ) = 

for some smooth function h{Z\,(x). Let z\ := y/\ai/~Re C| > so that z\ is the positive 
solution to the equation h(Zi,0) — 0. Then clearly Dz 1 h(zi,0) — — 2|ReC|zi 7^ and thus 
by the implicit function theorem there exists a smooth function Z\ : [0,/io) — > R>o of the 
form Zi(n) = z\ + O(fi) that solves h(Zi fi) — 0. Correspondingly, 

Si (A) = [Bi(A)| = Br(|A|5) : = Zt(\\\i) = y/\ ai X/RoC\ + 0(\\\) ~ y/\X\ = |A| K1 . 

is a branch of solutions to equation (|6.18|l . 

The imaginary part of equation (|6.17[) immediately gives that the frequency ui (A) is a 
smooth function of A, given by 

oj(X) = u + a 2 A + ImC*|Bi(A)| 2 + C(|A| 2 ) . (6.19) 

This finishes the analysis of equation (|6.16|l . 

Next, assume as induction hypothesis that we have found for 1 < j ' < n — 1, either for 
A G ( — Ao, 0] or for A G [0, Ao), solutions 

x (t) = 0,x 1 (t) = B 1 (X)e^ x)t ,...,x J (t) = B J (A)e^ (A)t 

of the equations ±0 = f(xo, • • • , %o', A), . . ., Xj = f(xj, . . . , xo; A) so that for all 1 < i < j, 
Bi(\) is a smooth function of |A| Ki , i.e. Bi(\) — _B*(|A| Ki ) ~ |A| Ki for some smooth function 
B* : [0, Aq ; ) — > C. We already proved this induction hypothesis for j = 1. 

Then by solving Xj+i = f(xj+i, . . . , xq\ A) we shall try to obtain Xj+x(t), that is we solve 

, 71—1 

=^F i (...;A)a; Inax0+1 _ ij o} = («^o + «A + C\x ]+1 \ 2 )x J+1 + p Xj + OQXp' 1 ) . 

i=0 

Now the Ansatz Xj+i(i) — Bj+ie 1 "'"* 1 '' leads to the equation 

iu(X)B j+ i = (ibJo+aX + ClBj+^Bj+i+pBj+OQX^- 1 ). (6.20) 

Anticipating that the solution Bj+i will have amplitude |A| K -' +1 , let us at this point define 

1 

the rescaled parameter /i := \X\ Kj+1 = |A| 2 3J and the rescaled unknown Bj+i =: (iZj+i. 
Then it holds that /j, 33 + = jAj' 1 ', which inspires us to define also the (smooth) rescaled 
functions Z* , . . . , Z* by 

Sr(|A| K1 ) =: M 3 ^i*(M 33 ), ■ • ••/-',:: A' . =: ^Z*^) for n = A^ 1 . 
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Then, because lj(X) = ujq + 0(A), equation (|6.20p obtains the form 



(c\z 3+1 \ 2 z 3+1 + pzM) + o(m 4 ) = o . 

Dividing by /i 3 we now find that we need to solve an equation of the form 

h(Z i+1 ,n) = C\Z J+1 \ 2 Z j+1 + pZi(Q) + 0{ji) . 

We solve this equation as follows. First of all, there is a unique Zj+i £ C for which 
h(zj+i,0) — C\zj +1 \ 2 Zj +1 + (3Zj(0) — 0. It clearly holds that z j+1 / because j3 / 
and Zj(Q) 7^ by the induction hypothesis. As a consequence, 

Dz j+1 h(z j+1 ,0) : v n- 2C\z j+ i\ 2 v + Cz 2 +1 v 

is invertible, since det Dz +1 h(zj+i, 0) = 3]C| 2 \zj+i 4 7^ 0. Thus, by the implicit function 
theorem, there exists a smooth function Z* +1 : [0, /10) — > C so that Z* +1 (0) = Zj+i 7^ and 
h(Zj +1 (n), fi) = 0. Correspondingly, 

B j+1 (X) = |B J+1 (A)| = B;+i(|A|5) := \\\ K ^z; +1 (\X\ K ' +1 ) ~ 



is a branch of solutions to (I6.20|l . This finishes the induction and proves the existence of the 
Hopf branch 

x (t) = 0, xi(t) = B!(A) e ^ (A)i , . . . ,x n (t) = B n (X)e^ Wt . 

The remaining branches in the statement of the theorem exist by symmetry. 

Finally, we consider the linearization of the normal form flow around the periodic solution 
on the r-th branch. Thus, we perturb our relative equilibrium by substituting into the normal 
form equations of motion 

x (t) = ey {t),. . . , x r -i(t) = ey r -i(t) 
x r (t) = (B 1 (X) + ey r (t))e^ (x)t ,...,x n (t) = (B„_ r+ i(A) + ey n (t))e^ X)t . 
This yields that y = M(X)y + 0{e), the linearization matrix being of lower triangular form 



/ 60(A) + ... + o»(A) 
* 60(A) 



M(A) = 







60(A) 
* bi(A) 



V 



bn— r+l(A) / 



The stability type of the Hopf curve is thus determined by the maps on the diagonal. The 
first of these, 60(A) + . . . + 6„(A) = ao(0) + 61 (0) + . . . + 6 n (0) + 0{\ X\), is hyperbolic because 
oo(0) + ... + fln(0) is hyperbolic and has the same eigenvalues as Oo(0) + 6i(0) . . . + a„(0). 

The maps So (A) = jojo + ctX + C(|A| 2 ) (r — 1 times) are hyperbolic for A / because 
on 7^ 0. Finally, the maps bi(A), . . . , b n _ r +i(A) are given asymptotically by 

bj{X)v = 2qs i | 2 w + CBjV + OflAr^M) . 

These maps are hyperbolic for A because detb^(A) = 3|C| 2 |Bj| 4 + OQXf****- 1 ) > 

□ 



and trbj(A) = 4(Re C) \B 3 \ 2 + 0(\X[ 



^ 0. Thus, the branches are hyperbolic. 



Remark 6.4 For any one of the Hopf branches given in Theorem 16.11 to be stable, it is 
necessary that ReC < 0. In turn, this implies that the branches exists for A with a\X > 0. 
For such A though, the eigenvalues of 60(A) = iujq + aX + C(|A| 2 ) have positive real part. 
Thus, the only branch of periodic solutions that can possibly be stable is the branch with 
the least synchrony and the largest amplitude, i.e. the branch with asymptotics 

x = 0,xi ~ |A| K1 , . . . ,x n ~ |A| K ™ . 

Indeed, this branch is stable precisely when ReC < and Uo(0) + ... + a„(0) only has 
eigenvalues with negative real parts. A 
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